Noise-resistant entanglement of strongly interacting spin systems 
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We propose and analyze a scheme that makes use of interactions between spins to protect certain corre- 
lated many-body states from decoherence. The method exploits the finite energy gap of properly designed 
Hamiltonians to generate a manifold insensitive to local noise fluctuations. We apply the scheme to achieve 
decoherence-resistant generation of many particle GHZ states and show that it can improve the sensitivity in 
precision spectroscopy with trapped ions. Finally we also show that cold atoms in optical lattices interacting via 
short range interactions can be utilized to engineer the required long range interactions for a robust generation 
of entangled states. 

PACS numbers: 



Quantum entanglement has recently emerged as an impor- 
tant resource in quantum information science and metrology 
Jll 01 . For example, entangled many -particle states can be 
used to perform long-distance quantum communication and 
scalable quantum computation and to enhance the spectro- 
scopic sensitivity in quantum-limited measurements. How- 
ever, many-particle entangled states are difficult to prepare 
and to maintain since they are extremely fragile: in practice, 
noise and decoherence rapidly collapse them into classical sta- 
tistical mixtures. 

In this Letter we propose and analyze a new method for the 
robust generation of entangled states and its protection against 
decoherence. Our approach is based on creating a degener- 
ate many-body subspace, which can be isolated from the rest 
of the Hilbert space by properly designed interactions. We 
focus specifically on interacting two-level systems where the 
two states are associated with spin sublevels, e.g trapped ions 
or neutral atoms yl IfU, |5l la \M ■ Such systems are of partic- 
ular importance for applications in precision measurements. 
Few-particle entangled states have already been prepared with 
trapped ions, and proof-of-principle experiments demonstrat- 
ing the improvement of spectroscopic sensitivity have been 
carried out QHl- Our aim is to devise a technique to increase 
the robustness of entangled collective states. We discuss po- 
tential applications of this protected evolution to perform pre- 
cision measurements with trapped ions and to generate many- 
particle-GHZ-type states in optical lattices. 

Our approach can be best understood by considering a sim- 
ple example a multi-spin system with isotropic ferromagnetic 
interactions. These interactions will naturally align the spins. 
While all of the spins can be rotated together around an ar- 
bitrary axis without cost of energy, local spin flips are ener- 
getically forbidden. Consequently ferromagnetic interactions 
enable the generation of superpositions (suitable for exam- 
ple for precision spectroscopy) with a substantial suppression 
of decoherence (See Fig.l). We note that in Ref.|8|] gapped 
quantum systems have been proposed for decoherence-free 
entanglement generation by adiabatic ground-state transitions. 
Specifically, Ref. |8[] exploited the finite energy gap to protect 
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FIG. 1: Many-body protected manifold (MPM) generated by the 
isotropic ferromagnetic exchange: (a) The ground manifold is 
spanned by the states with all spins aligned in the same direction, (b) 
Any global rotation about an arbitrary axis is allowed, (c) local spin 
flips are energetically suppressed. Panel (d) shows that the MPM lies 
on the Bloch surface with radius J = Jm ax — N/2. 

non-degenerate ground states. The present mechanism is ca- 
pable of protecting degenerate multi-state manifolds with a 
large number of quantum degrees of freedom and use non- 
equilibrium dynamics as the mechanism for entanglement 
generation. 

We consider a collection of TV spin 1/2 atoms with isotropic 
infinite range ferromagnetic interactions described by the 
Hamiltonian H = H pro t + H m , where 

H P rot = -\J m and H m =w j| 0) . (1) 

We used Ja to denote collective spin operators of the N 
atoms: Ji "* = h af, where a = x, y, z and af is a Pauli 
operator acting on the i th atom. Moreover we have identified 
the two relevant internal states of the atoms with the effective 
spin index a = Til- These states have energy splitting ujq 
in units such that K = 1. 

An appropriate basis to describe the dynamics of the system 
is the basis spanned by collective pseudo-spin states denoted 
as | J, M, [3} z |0]. These states satisfy the eigenvalue relations 
J(°) 2 | J, M, 0) z = J(J + 1)| J, M, p) x and J { Z 0) \J, M, f3) z = 
M\J, M, (3) x , with J = N/2, . . . ,0 and —J < M < J. f3 
is an additional quantum number associated with the permu- 
tation group which is required to form a complete set of labels 
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for all the 2 N possible states. 

The isotropic Hamiltonian H prot has a ground state man- 
ifold spanned by a set of TV + 1 degenerate states. They 
lie on the surface of the Bloch sphere with maximal radius 
J = N/2 (see Fig.l) and are totally symmetric, i.e. invari- 
ant with respect to particle permutations. This set of states 
are fully characterized by J and M, and we denote them as 
| N/2, M) z with the additional label j3 omitted. There is a fi- 
nite energy gap E g = XN that isolates the ground state man- 
ifold from the rest of the Hilbert space. This gap is the key 
for the many-body protection against decoherence. We will 
refer to the ground state manifold as the many-body protected 
manifold (MPM). 

To understand the protection within the MPM, we first as- 
sume that the dominant type of decoherence is single-particle 
dephasing. Such dephasing comes from processes that, while 
preserving the populations in the atomic levels, randomly 
change the phases leading to a decay of the off-diagonal den- 
sity matrix elements. We model the phase decoherence by 
adding to Eq. (1) the following Hamiltonian ifioll 



H env = \^h i {t)a z i 



(2) 



where the hi(t) are assumed to be independent stochastic 
Gaussian processes with zero mean and with autocorrelation 
function hi(t)hj(r) = 6ijf(t — r). Here the bar denotes av- 
eraging over the different random outcomes. In what follows 
we will use the property that zero mean Gaussian variables 

satisfy the property exp[— i J Q drh(T)] = exp[—T(t)], with 

Let us first study the A = dynamics, with the sys- 
tem in the ground manifold at t = 0. Phase decoherence 
modifies the free precession of the Bloch vector generated 
by H m . While (Jz (i)} remains a conserved quantity, the 
x and y projections decay exponentially with rate T(t), i.e. 

&)) = e- r M(ji°](t))| r =o, where (ji 0) (i)>|r=o is the 
expectation value in the absence of noise ifioll . In addition, 
due to local phase fluctuation which drive the system out of 
the ground state manifold and deplete the J = N/2 levels 
also ( j( ' 2 ) decays exponentially. 

The effect of the environment on the free evolution is sig- 
nificantly reduced by H prot . The protection can be best un- 
derstood by using the basis of collective states. In terms of 
collective spin operators H env can be written as: 
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where g k {t) = -^^ =1 h 3 (t)e 
5 2_(5=i a 3 e% N ■ I n tne presence of a large energy gap E g 
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one can distinguish two different types of processes: (i) Deco- 
herence effects that take place within the MPM due to the col- 
lective dynamics induced by the k = component of H env , 
and (ii) transitions across the gap induced by the inhomoge- 
neous terms with k ^ 0. The latter couple the MPM with the 



rest of the system. All allowed transitions must conserve M 
since both the system and noise Hamiltonian commute with 

f(o) 

J z 

In the limit when the noise is sufficiently slow ( i.e. when 
the spectral density of the noise has a cutoff frequency ui c <C 
Eg) the effect of phase decoherence is dramatically reduced 
as type (ii) processes are energetically forbidden and only 
type (i) processes are effective. In this limit the noise acts 
as a uniform random magnetic field. If at t = the density 

matrix, p, is p = J2m,m />m,m(°) If > M )(f > then af " 
ter time t each component p M A ~ 7 acquires an additional phase 
e i{9 M {t)-9 a {t)) with q m (q = 7= Jo" 9°(t), such that 



Pm,mW = PM.Al(°) e 



mj (M-M)t _-r(t)i 



(4) 

Note the factor of s/N in the denominator of 9m- It 
is fundamental for the reduction of the effect of decoher- 
ence within the MPM. For example, Ji°l decays N times 
slower than in the unprotected system: i.e. (Jxy(t)} = 

e- r ^(J^(t))|r=o. 

Let us now elaborate more on the necessary conditions re- 
quired to suppress type (ii) processes. For this, we will assume 
the power spectrum of the noise, f(uj) = J dte~ tut fit), to 
have a cut-off frequency uj c (e.g. /(w) = / for uj < uj c 
and otherwise with / < A ). A time dependent pertur- 
bation analysis predicts that as long as oj c < E g , the decay 
rate j M (t) of the diagonal elements pM iuif) of the den- 
sity matrix due to type (ii) process is always bounded by 

"*w < (^) (#0 



7 



<C 1. Hence in the slow noise 



limit the population of the MPM is fully preserved. 

Protected generation of N-particle GHZ states: A^-particle 
Greenberger-Horne-Zeilinger (GHZ) states ll ill of the form 

N _N - 
1 2 ' 2 } 



GHZ 



~ V2\ l 2'2 )x + C 



where \ ^j,±^-) x are fully polarized states along the ±x and 
<f>± are arbitrary real phases, can be generated by evolving the 
state \N/2, N/2) x with an interaction of the form ILlfllll 

Hz= X J { z 0)2 - (5) 
The generation process can be understood by writing 
\N/2,N/2) X as Y jM C m \N/2,M) z . During the evolution 
the Hamiltonian imprints an M 2 dependent phase to the dif- 
ferent components. As the system evolves, at first the winding 
of the phases leads to a collapse of (Jx }■ However, at time 
X^rev = 1" all components rephase with opposite polariza- 
tion, and a perfect revival of the initial state is observed with 
(jj, 0) ) = -N/2 (see Fig. 2). Right at time t = t rev /2 the 
system becomes an Af-particle GHZ state. 

Recent experiments 0] have used this type of scheme 
to generate GHZ states in trapped ions with the aim to per- 
form precision measurements of uiq ( ideally the use of GHZ 
states should enhance the phase sensitivity to the fundamental 
Heisenberg limit IU4lD . However, decoherence significantly 
limited the applicability of the method and in practice even 
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for only six ions, the achieved phase sensitivity was signifi- 
cantly below the fundamental limit. The effect of decoherence 
can be quantified by calculating the fidelity of creating a GHZ 
state defined as T(t ) = (ip§ HZ \p(t Q )\^ HZ ). It is possible 
to show that, after the evolution with the total Hamiltonian, 
Henv + H z , the fidelity is degraded to: 

/ 1 + e -r(t )\ N 
^(*o)= ^ ■ (6) 



The TV-dependent decay of J 7 (to) and the exponential de- 
cay of ( Jx ) reflect the rapid collapse of the entanglement as 
TV increases (See Fig. 2). Consider now the same generation 
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FIG. 2: (Color online) Fidelity of GHZ generation vs. TV with and 
without gap protection. We assume a noise with step-like spectral 
density with amplitude / and cutoff frequency u) c = X ■ The dashed 
green and solid red lines are for a protected system, A = 50, with 
/ = 0.6x an d / = O.lx respectively. The dot-dashed blue and 
black dotted lines are for an unprotected system, A = 0, and same 
decoherence parameters. In the inset we show (Jx°\t)) /TV for TV = 
50. 

process within the MPM. The latter can be realized by replac- 
ing H z — > Hprot + H z . In the slow noise limit the reduced 
density matrix evolves as Eq. (0]i with the phase ujo(M — M) 
replaced by x(M 2 — TV/ 2 ). For the initially polarized state the 
fidelity is given by: 

The independence of J- (to) on TV, and the N-times reduced 
decay rate of (Jx ) demonstrate the usefulness of MPM to 
generate a large number of entangled particles. 

Let us now discuss physical implementations of the MPM. 



In Ref. [81 it has been shown that the Hamiltonian H, 



prof ' 



can be implemented by using the collective vibrational mo- 
tion of the ions in a linear trap driven by illuminating them 
with a laser field. If the detuning of the laser, 6, from the in- 
ternal transition frequency is large compared to the linewidth 
of the resonance but sufficiently different from the ion vi- 
brational frequency, then the dominant processes are two- 



photon transitions which lead to simultaneous excitations of 
pairs of ions and thus to an effective Hamiltonian of the form: 
H ef f = x(J m ~ Ji° )2 ) + (3Jz 0) - The cost of implementing 
H e f f instead of H z is the additional echo techniques required 
for making the former insensitive to heating of the vibrational 
motion of the ions: while H z is naturally insensitive to ther- 
mal vibrations due to the fact that it is generated by bichro- 
matic beams lfl2h ( which eliminate any dependence on vibra- 
tional quantum numbers due to the destructive interference be- 
tween the transition paths) H e f f lacks this insensitiveness as 
it is generated by monochromatic beams. On the other hand, 
the implementation of H e j f has the advantage in comparison 
with H z that the MPM protects the system from local pertur- 
bations caused by non-ideal conditions, such as differences in 
the Rabi frequency experienced by the string of ions, or by ex- 
ternal disturbances, such as magnetic field noise. It should be 
noted, however, that the system is not protected by the MPM 
against spontaneous emission and global perturbations. 

MPM with local interactions: Up to now we have explored 
only the generation of an MPM via isotropic long-range in- 
teractions. In practice, however, it is desirable to have a sim- 
ilar kind of protection generated by systems with short range 
interactions such as those provided by cold atoms in optical 
lattices. We now show how an MPM can be created in lattice 
systems and can be used to robustly generate TV -particle GHZ 
states in these systems. We consider ultracold bosonic atoms 
with two relevant internal states confined in an optical lattice 
with unity filling deep in the Mott insulator regime. Such sys- 
tems is described by an effective spin XXZ Hamiltonian 11511 : 



H lat = Hh + Hi = -A E °i & J 



<i,j>,a 



x E w w 

<i,j> 



Here we have identified the two possible states at each site 
with the effective spin index a =|, J,. The sum of (i,j) is over 
nearest neighbors. In Eq. dH) the coefficients are A = r 2 /U-[i 
and x = t2 (U\i + — 2U-\i)/U 2 , where U is the mean 
of the three different spin dependent on-site interactions ener- 
gies JTfr, Un, U-\i, which we assumed only slightly different 
between each other, and r is the tunneling energy between ad- 
jacent sites which we assume spin independent. Both U and 
r are functions of the lattice depth. For simplicity we restrict 
the analysis to one dimensional systems and assume periodic 
boundary conditions. 

Hh is spherically symmetric and in terms of collective spin 
operators it can be written as 

k>0,a v 7 

All the TV + 1 states with J = TV/2 are degenerate and span 
the ground state of Hh thus form an MPM. Hi is not spheri- 
cally symmetric but we can also write it in terms of collective 
operators as 



Hi = -%n o)2 -Se j^j^ cos ( ^ 

TV TV ^ V TV 

fe=i v 



(10) 
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FIG. 3: (color online) a) Fidelity to generate a GHZ state vs X/x- 
In the inset we show ( J^ ' {£)}. The blue dot-dashed, dotted black, 
dashed green and solid red correspond to A = 0, 5, 10, 20 respec- 
tively. The plots are obtained by numerical evolution of Eq. {8} for 
N = 10. b) In the presence of phase decoherence the fidelity of the 
GHZ state preparation in the lattice is degraded as iV grows because 
in the lattice the gap decreases and the generation time increases with 
N. In this plot we used a system with lo c — x, A = 100x and 
T = O.Olx- At N = 90, AE g = uj c and it explains the drop of T 
for N > 90. 

If the condition x -C A is satisfied, which is naturally the case 
for spin independent lattices, the effect of the Ising term can 
be treated perturbatively. Assuming that at t — the initial 
state is prepared within the J = N/2 manifold, a perturbative 
analysis predicts that for times t such that yt < X/x, Hh 
confines the dynamics to the MPM and transitions outside can 
be neglected. As a consequence, only the projection of Hi 
on Hh is effective. As VHi = XeJz°^ 2 + const, with Xe = 



■j^filM, Hj acts as a long range Hamiltonian. 

In Fig. 3a we show the fidelity to create a GHZ state at 
Xet = 7r/2vs X/x and contrast the dynamical evolution of 
(Jx^) for different X/x ratios assuming at time t — all the 
spins are polarized in the x direction. If only the Ising term is 
present, A = 0, it induces local phase fluctuations that leads 
to fast oscillations in (J| 0) ) = N/2 cos 2 [2xt]. On the other 



hand, as the ratio X/x increases, the isotropic interaction in- 
hibits the fast oscillatory dynamics and instead ( Jx ) exhibits 
slow collapses and revivals. For A 3> \ me dynamics exactly 
resembles the one induced by H z and at Xet = 7r/2 the initial 
coherent state is squeezed into a GHZ state. 

Hff also provides protection against phase decoherence. 
However, Hh is not as effective as H prot because the energy 
gap between the MPM and the excited states of Hh (EqfTOl) 
vanishes in the thermodynamic limit as E g — ► X/N 2 . This 
is a drawback of the short range Hamiltonian for the purpose 
of fully protecting the ground states from long wave length 
excitations [19]. Nevertheless, the many body interactions 
can still eliminate short-wavelength excitations since in the 
large N limit they remain separated by a finite energy gap, 
8A. In Fig. 3b we quantify the protection provided by the lat- 
tice Hamiltonian against population decay by plotting J 7 as a 
function of N . As expected, an abrupt drop of the fidelity 
occurs at the value of N at which E g = uj c . 

In summary we proposed and described a new method for 
the robust generation of entangled states and for their protec- 
tion against decoherence and discussed its applicability for the 
generation of many-particle-GHZ-type states in ion traps and 
optical lattices. We emphasize that, even though we have lim- 
ited the discussion to ensembles of spin S = 1/2 particles, the 
MPM ideas can be straightforwardly generalized to systems 
composed of higher spin atoms. Besides entanglement gener- 
ation, the MPM might have also important applications for the 
implementation of good storage memories using for example 
nuclear spin ensembles in solid state 

El 

or photons jl7ll . 
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